Harmonic oscillator chains as Wigner Quantum Systems: 
periodic and fixed wall boundary conditions in g[(l|n) solutions. 



We describe a quantum system consisting of a one-dimensional linear chain of n identical 
harmonic oscillators coupled by a nearest neighbor interaction. Two boundary conditions are 
taken into account: periodic boundary conditions (where the nth oscillator is coupled back 
to the first oscillator) and fixed wall boundary conditions (where the first oscillator and the 
nth oscillator are coupled to a fixed wall). The two systems are characterized by their Hamil- 
tonian. For their quantization, we treat these systems as Wigner Quantum Systems (WQS), 
allowing more solutions than just the canonical quantization solution. In this WQS approach, 
one is led to certain algebraic relations for operators (which are linear combinations of position 
and momentum operators) that should satisfy triple relations involving commutators and anti- 
commutators. These triple relations have a solution in terms of the Lie superalgebra g[(l|n). 
We study a particular class of g[(l|n) representations V{p), the so-called ladder representations. 
For these representations, we determine the spectrum of the Hamiltonian and of the position 
operators (for both types of boundary conditions). Furthermore, we compute the eigenvectors 
of the position operators in terms of stationary states. This leads to explicit expressions for 
position probabilities of the n oscillators in the chain. An analysis of the plots of such position 
probability distributions gives rise to some interesting observations. In particular, the physical 
behavior of the system as a WQS is very much in agreement with what one would expect from 
the classical case, except that all physical quantities (energy, position and momentum of each 
oscillator) have a finite spectrum. 
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1 Introduction 



Coupled systems describing interaction of oscillating or scattering subsystems and the correspond- 
ing operators have been widely used in classical and quantum mechanics. For example, coupled 
harmonic oscillators have been studied in quantum information theory, quantum optics (photonic 
crystals), and for describing phonons in a crystal [1-6]. 

In a previous paper [7] , we initiated the analysis of such systems as a Wigner Quantum System 
(WQS) [8-10]. The system studied in [7] consists of a string or chain of n identical harmonic 
oscillators, each having the same mass m and frequency to. The position and momentum operator 
for the rth oscillator (r = 1,2, ... ,n) arc given by Qj- and pr] more precisely Qj- measures the 
displacement of the rth mass point with respect to its equilibrium position. The oscillators are 
coupled by some nearest neighbor coupling, represented by terms of the form (g^ — Qr+i)'^ in the 
Hamiltonian. In [7], we considered such a system with periodic boundary conditions, where the 
last oscillator is again coupled to the first one (i.e. Qn+i = Qi)- In other words, the oscillators are 
on a circle rather than on a line. The Hamiltonian of this system will be denoted by Hp. In the 
present paper we will reconsider this system, and at the same time also concentrate on a second 
type of boundary conditions: fixed wall boundary conditions. For this second case, the first and 
last oscillator are coupled to a fixed wall (i.e. = Qn+i = 0); the Hamiltonian of this second system 
will be denoted by Hpw- 

The original contribution of [7] consisted of the treatment of the system described by Hp as a 
Wigner Quantum System [8-10]. This quantization procedure is based upon the compatibility of 
Hamilton's equations and the Heiscnberg equations. In this procedure, one does not require the 
canonical commutation relations (CCRs), but instead uses more general relations describing the 
before mentioned compatibility. As a consequence, one obtains different classes of solutions for 
the system, of which the canonical solution is only one. For a more detailed description of WQSs, 
see [7] and references therein. The main feature of this treatment is that the CCRs are replaced 
by these compatibility conditions (CCs). For the Hamiltonian Hp, these CCs are triple relations 
in terms of operators (r = 1,2, ... ,n), involving both commutators and anti-commutators. 
Herein, the operators af" are certain linear combinations of the position and momentum operators. 
In [7] we showed that these triple relations have a solution in terms of certain generators of the Lie 
superalgebra fl[(l|n) [11]. Furthermore, we considered a class of Fock representations [12] of gl(l|n) 
and analysed the energy spectrum and the position operator spectrum in these Fock representations 
W{p). These Fock representations are rather restricted, however, and do not illustrate the features 
of general unitary irreducible representations [13] of g[(l|n) (more precisely, of its compact form 
u(l|n)). 

In the present paper, we consider a more general class of st(l|n) representations, the so-called 
ladder representations V{p) [13]. These representations are also easy to describe, but more im- 
portantly they show interesting properties of the physical quantities (energy spectrum, position 
operator spectrum) of the system, far more general than those of the Fock representations W{p). 

A second original contribution of this paper is that we show how also the second system (fixed 
wall boundary conditions) described by Hpw can be treated as a Wigner Quantum System. In 
fact, we show that at the algebraic level the two Hamiltonians give rise to the same triple relations, 
but with different constants. As a consequence, the two systems can be treated similarly, both 
having solutions in terms of the Lie superalgebra fl[(l|n). 

The contents of the paper is as follows. In Section 2 we shortly review the treatment of Hp as 
a WQS [7], and then analyse the Hamiltonian Hpw of the fixed wall boundary case in a similar 
way. As the algebraic triple relations are, up to different constants, the same, we describe the 
gl(l|n) solution in a subsection. Special attention is paid to determining the critical value for the 
coupling constant c, for which a 0t(l|n) solution exists in the fixed wall boundary case. We also 
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describe briefly the new class of representations V{p), the ladder representations, and the explicit 
action of 3[(l|n) generators on simple basis vectors w{9;s) of V{p). In Section 3 we determine the 
spectrum of the Hamiltonians Hp and Hpw in the representations V{p). Due to the fact that these 
Hamiltonians have a diagonal action in the current model, their spectrum is easy to determine. We 
compare the two cases, and discuss some aspects of degeneracy of the energy levels. 

The mathematically more difficult problem is that of determining the eigenvalues and eigenvec- 
tors of the position operators qr, for these representations V{p), since qr is a general odd operator 
of 0l(l|n) with a nondiagonal action on the standard basis of V{p). Fortunately, we have managed 
to give a general treatment (and procedure) of how to construct such eigenvalue spectrum and 
eigenvectors for an arbitrary odd element of 0l(l|n) in an arbitrary unitary representation in [14]. 
Using the techniques of [14], we manage to determine the eigenvalues of qr for both systems in 
Section 4, and to construct the eigenvectors of qr in Section 5. In fact, the eigenvalues of qr for the 
system described by Hp were already given in [14] as an example of the new technique, so only 
those for the second system described by Hpw are new. We also give some plots of the spectrum 
of the position operators, as a function of the coupling constant c, and make some observations 
related to these plots. 

Section 5 is devoted to computing position probability distributions for the two systems under 
consideration. Since the techniques of [14] allow us to give an explicit expansion of the normalized 
eigenvectors of qr, for any of its eigenvalues ix^, in terms of the stationary states w{9;s), we can 
invert these relations and express the stationary states w(9;s) in terms of the eigenvectors of ijr 
(for any particular r). The square moduli of the coefficients in these expressions have the usual 
quantum theory interpretation as the probability of finding the rth oscillator in position 
when the system is in its stationary state w{9;s). We compute these probabilities analytically 
(for small p- values), and numerically for some examples. We also give a number of plots of such 
probability distribution functions. These yield the "spacial properties" of the oscillators. The 
results are quite interesting, and show some analogy of what one would expect of the system from 
a classical treatment. For the ground state (stationary state of lowest energy), the position with 
highest probability (for each oscillator) is the equilibrium position (corresponding to Xp = 0), with 
strongly decreasing probability as the (discrete) eigenvalue is further away from 0. For the most 
excited state (stationary state of highest energy), the position with highest probability is away 
from (symmetrically to the left and the right). The higher the coupling constant, the further 
away this highest probability position. The probability of the equilibrium position is zero here. All 
oscillators have the same position probability distribution in the periodic boundary CcLSG, cLS they 
are all equivalent (completely symmetric on a circle). In the case of fixed wall boundary conditions, 
the behavior also follows the classical properties, with the first and last oscillators (those fixed to 
the wall) "oscillating less" compared to those away from the wall. We end the paper by some 
concluding remarks in Section 6. 

2 The Hamiltonian of the system: solutions as a WQS 

2.1 System with periodic boundary conditions 

In this subsection, we briefiy recall the algebraic treatment of the Hamiltonian for a system con- 
sisting of coupled harmonic oscillators, with periodic boundary conditions, treated as a Wigner 
Quantum System [7]. In several models [2-6] such a quantum system consisting of a linear chain 
of n identical harmonic oscillators coupled by springs is used. The Hamiltonian of such a system 
is given by: 

-2 2 

Hp=j:[t^ + ^^^^ - )' ^^-'^ 

r=l 
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where each osciUator has mass m and frequency oj, q^. and p,. stand for the position and momentum 
operator for the rth oscillator (or rather, qr measures the displacement of the rth mass point with 
respect to its equilibrium position), and c > is the coupling strength. In the case of periodic 
boundary conditions (indicated by the index in Hp), one assumes in (12. ip 



q-n+i = ^1- (2.2) 

With these periodic boundary conditions one can think of the oscillators as being located on a 
circle, as the last oscillator is again coupled to the first one. 

In [7], it was shown that one can relax the canonical commutation relations for the operators 
qj. and p^, leading to a larger class of solutions for the system described by (|2.1|) . This is known 
as a Wigner Quantum System approach [9, 10]. In this approach, one imposes the compatibility of 
Hamilton's equations 

, dHp dHp 

qr = -^^ Pr = — WT- [r = l,2,...,n) (2.3) 

OPr Oqr 

(formal derivatives) and the Heisenberg equations 

Pr = -^[Hp,pr], qr = -^[Hp,qr] (r = 1, 2, . . . , n) (2.4) 
when viewed as operator equations. These compatibility conditions (CCs) read 

[Hp,qr] = Pr, (2.5) 

m 

[Hp,Pr] = —ihcm (jr^i + ihm{uj'^ + 2c) qr — ihcmqr+i, (2-6) 

with r £ {1,2,... ,n}, (jn+i = qi and % = qn- Under the canonical commutation relations, the 
CCs (|2.5p and (|2.6p are automatically satisfied. But the system (j2.5p - (|2.6p . with Hp given by (|2.ip . 
has other interesting solutions [7]. 

To study these other solutions, one introduces discrete Fourier transforms of the (self-adjoint) 
operators and pr by 




^ ' ^-2nijr/n^+ ^ ^2nijr/n- | ^ ^2.7) 



2mnujj \ ^ ^ 



Pr = yj\j (e-2-i-/-a+ - e2-^'^/"a- ) , (2. 



where ujj are positive numbers with 



oj] = + 2c - 2ccos(^) = + 4csin2(^), (2.c 



n n 



and are operators satisfying (a^)^ ~ ' terms of these new operators, the Hamiltonian 
reads [7] 

^p = E^(«74 + 4«7)- (2-10) 

Note that we no longer require the canonical commutation relations for the operators qr and pr, so 
also the operators no longer satisfy the usual boson relations [a^,a^] = and [aj,a^] = 5jk- 
In the WQS approach, the relations that should be satisfied follow from (j2.5p - ()2.6p . and read 
explicitly [7]: 

n 

\^u;j{a^a+ + a+aj),a^ =±2uJka^, (fc = 1, 2, . . . , n). (2.11) 
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2.2 System with fixed wall boundary conditions 

In this subsection we consider a similar system, again consisting of coupled harmonic oscillators, 
but this time with fixed wall boundary conditions. Now the Hamiltonian reads: 



HfW = E (£: + -^"Ir + -irilr - qr+lf), (2.12) 



V2m ■ 2 '''' ' 2 

r=l 

with the same data as in ()2.ip . but 

^0 = qn+i = (and po = Pn+i = 0). (2-13) 

In other words, we assume that the first and last oscillator (i.e. the oscillators numbered 1 and n) 
are attached to a fixed wall. 

The treatment of this Hamiltonian in the WQS approach is very similar to the previous sub- 
section. In fact, the CCs (j2.5p - (|2.6p remain the same (with Hp replaced by Hfw^ but with 
Qo = Qn+i = 0). These different boundary conditions lead to a different type of transform. Instead 
of a discrete Fourier transform, we now use a discrete sine transform. More explicitly, we introduce 
the following transformations of the (self-adjoint) operators qr and pr'- 



where the uij are positive numbers given by 



= t^2_^2c-2ccos(^^) = u;^ -h4csin^( , J . (2.16) 



The operators a- satisfy the adjointness conditions 



In terms of these new operators, the Hamiltonian ()2.12p is given by: 



(a±)t = a^. (2.17) 



^F"^ = E ^ ("7 4 + ) ' ^2.18) 



and the compatibility conditions become 

n 

[j^cD,(a-a+ + a+aT),a±] =±2wfca± (/c = 1, 2, . . . , n). (2.19) 

i=i 

So the algebraic expression of the Hamiltonian, and the CCs, are the same in the two cases 
considered, apart from the replacement lOj — > tDj. This implies that the algebraic solutions will be 
similar, even though the conclusions about physical properties will be different due to the different 
numerical values of the numbers ujj and iby 
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2.3 The solution 

In was shown in [7] that the triple relations (j2.1ip involving both anti-commutators and commu- 
tators have a solution in terms of generators of the Lie superalgebra 0[(l|n) [11]. More explicitly, 
let £|l(l|n) be the Lie superalgebra with standard basis elements ejk {j, k = 0,1, . . . ,n) where e^o 
and eofc {k = 1, . . . ,n) are odd elements (degeo/t = degCfco = 1) and the remaining basis elements 
are even (having degree 0), with bracket 

N, efc^l = 6,ke^l - (-l)d<^s(-».)deg(e.05^,e, (2.20) 
and star condition ej^- = eji. Then a solution of (j2.1ip is provided by 




a-j = ejo, aj = eoj {j = l,...,n) (2.21) 



where 

1 " 

Pj = -^j + 7^'^'^' (J = l,---,^)- (2-22) 

n — 1 ^-^ 

k=\ 

All these numbers /3j should be nonnegative. By the periodic boundary conditions, the /3j's satisfy 
fin-j = Pj, and 

Pi > P2> ■■■ > P\n/2\ 7 /5[n/2j < P\n/2\+l < ' ' ' < Pn- (2.23) 

It was analysed in [7] that all these /3j's are indeed positive provided the coupling constant c lies 
in a certain interval [0, co[, with cq some critical value depending upon n. 

In the case of fixed wall boundary conditions, the analysis is slightly different. A solution 
of p.igp is given by 





aJ = \^ejo, 4 = \/^eoi {j = h---,n) (2.24) 



where in this case 



Pi = -u 



1 " 

'j + 7y2^k, (i = l,...,n). (2.25) 

n — 1 ^-^ 

k=l 

Again, all these numbers Pj should be nonnegative. First, note that for c > one has that 
Oi < UJ2 < ■ ■ ■ < i^n and hence that 

Pi> p2> ■■■> Pn. (2.26) 

Thus all Pj are positive if and only if /?„ is positive. Secondly, for c = one has that Pj = 
Lo/{n — 1) > 0, and since Pn is a continuous function of c there exist positive values of c such that 
Pn > 0. Thus, there exists (in general) a critical value cq such that each Pj > for c < cq and such 
that for c = Co one has that /?„ = 0. The same upper bound on the critical value cq applies as on 
Co in periodic boundary conditions case since one can mimic the proof of [7, Proposition 2]. Since 
for 77. = 2 we have that P2 = + c there are no conditions on c in this case. Also for n = 3 
there are no conditions since one can verify (numerically) that in this case P3 > 0. In Tabled! we 
compare the critical values cq/w^ of the periodic boundary conditions case with the critical values 
cq/uj'^. One notices that these critical values are interleaved, and that they become more and more 
similar as the number of oscillators n increases. 
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n 






n 




co/uj'^ 


4 


0.9873724357 


2.1108888881 


13 


0.10546881460 


0.10521909714 


5 


0.7500000000 


0.7016444817 


14 


0.09256321610 


0.09509684206 


6 


0.3457442295 


0.4138598334 


15 


0.08687882025 


0.08675785013 


7 


0.2982653656 


0.2921798279 


16 


0.07814800074 


0.07976866442 


8 


0.2061705212 


0.2254893243 


17 


0.07388896853 


0.07382573538 


9 


0.1851128402 


0.1835156565 


18 


0.06760983697 


0.06871018095 


10 


0.1464642846 


0.1547079900 


19 


0.06429500840 


0.06426020587 


11 


0.1343028683 


0.1337254495 


20 


0.05957194222 


0.06035363583 


12 


0.1134651313 


0.1177656002 


21 


0.05691629341 


0.05689649085 



Table 1: Critical values cq/lo'^ (periodic boundary conditions) and cq/uj'^ (fixed wall boundary 
conditions) . 



2.4 A class of g[(l|n) representations V{p) 

In the case of canonical commutations relations, there is essentially only one representation of the 
system, following from the Heisenberg-Weyl algebra satisfied by the operators (jr and pr- In the case 
of WQS, the properties of the systems described by the Hamiltonians Hp and Hpw depend on the 
£|[(l|n) representation considered. In principle, any unitary representation of g[(l|n) can be taken 
into account [13]. A simple class of Fock representations W{p) [12] was already investigated in [7], 
for the case of periodic boundary conditions. This class of Fock representations is easy to work 
with, but also rather restricted as the basis vectors involve "fermionic" variables only, see [7, (4.1)]. 

In this paper, we will consider a richer class of representations, the so-called ladder repre- 
sentations V{p) [13]. These representations were considered as a special case in [14]. They are 
characterized by a positive integer p, and are finite-dimensional unitary representations atypical of 
type 2 [14]. A simple notation for the vectors of V{p) is: 

w{e; s) = w{9; si, S2, • • • , Sn), G {0, 1}, Si G {0, 1,2,.. .}, and + si + • • • + s„ = p. (2.27) 

Thus here the basis vectors involve one "fermionic" variable 9 and n "bosonic" variables Sj. In this 
notation the highest weight vector is w{l;p — 1,0, . . . ,0). 

The action of the 0[(l|n) generators on the basis ()2.27|) is given by (1 < A; < n) [13, 14]: 

eoow{e;s) = 9 w{9;s), (2.28) 

ekkw{9; s) = Sk w{9; s), (2.29) 

ekQw{9; s) = 9^/sk + l w{l - 6*; si, . . . , + 1, . . . , s„), (2.30) 

eQkw{9; s) = (1 - 6')^/i^ w{l - 9; si, . . . , Sk - 1, ■ ■ ■ , s„). (2.31) 

From these one deduces the action of other elements Cki- The basis w{9;s) ol V{p) is orthogonal, 
i.e. {'w{9; s),w{9'; s')) = 5Q^gi5s^s', and with respect to this inner product the action of the generating 
elements satisfies the conjugacy relations e|.Q = eok and ej,^ = e^o- 

3 On the spectrum of the Hamiltonians in the ladder representa- 
tion 

In this section, we study the spectrum of the Hamiltonians (j2.ip and (j2.12p in the ladder represen- 
tation V{p) of the g[(l|n) solution. Although the basis vectors of the representation are in both 
cases eigenvectors of the Hamiltonian, the spectrum in the two cases is quite different. 
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3.1 Energy eigenvalues in case of periodic boundary conditions 

For the periodic boundary conditions case, the Hamiltonian ()2.10p is given, using (j2.2ip . by: 

n 

Hp = h{(3 epo + y^/3fc efcfc), 

k=l 

with (3 = Ylk=i f^k = Z]fc=i ^k- Since the action of each Skk is diagonal in the basis w{9; s), see ()2.28p 
and (j2.29p . this imphes that each basis vector w{6; s) is an eigenvector of Hp, or a stationary state. 
Indeed, one finds that 

n 

Hp w{9; s) = h{l3e + l3kSk)w{e; s) = hEe,s wiO; s). 

k=l 

When c = 0, one has that (3^ = uj/{n — 1) and (3 = Lon/{n — 1), so in this case there are only 
two eigenvalues namely 



with multiplicities 



Eo,s = 7 and Eis = + w 

n — 1 n — 1 



p + n—l\ fp + n — 2 

and 

n — 1/ \ n — 1 



respectively. 

When c > 0, these two levels each split into a number of energy levels with lower degeneracies. 
Recall that (5k = Pn-k- So, when n = 2r is even, one can rewrite Eq g as follows: 



r-l 



{Sk + Sn-k) + PrSr + PnSn- (3.1) 



k=l 

It is then clear that any basis vector w{6; s') for which s'^ + s^_^ = + Sn-k {k = 1, . . . ,n — 1) 
yields the same Hp eigenvalue, independent of the value of c, i.e. for such 6 and s' one has that 
Eg s = Eq s'- Alternatively, one can say that for a fixed value of 6, Eq g is completely determined 

by' 

(Sl +S„_l,S2 + ■571—2) • • • ) Sr— 1 + s„_,.+i,Sr,s„), (3.2) 

which is a composition of p — into r + 1 parts. In the same way one sees that when n = 2r + 1 is 
odd, Eg^g with 9 fixed, is determined by 

(Sl + S„-1,S2 + 

For < c < Co and n = 2r or n = 2r + 1, the number of different energy levels is thus in general 
given by 

P + ^\^(P + ^-^)_ (3.4) 



p J V p ~ 1 

One can also say something about the degeneracy of an individual energy level hEg^g. The 

number of compositions of an integer in two parts is N + 1. Using this fact and (|3.2p or (j3.3p 
one sees that the degeneracy of hEg^g is at least: 

L(n-1)/2J 

Yl (sk + Sn-k + l)- (3.5) 
k=l 
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Figure [J shows two spectra of Hp for n = 4 and n = 5 with p = 2 va. both cases. These figures 
confirm our findings: in general for < c < cq there are 

/2 + 2\ [2 + 2- l\ 

different energy levels. For n = 4, the energy level hEg^g with threefold degeneracy is given by 

-E'0;2,0,0,0 = -E'0;0,0,2,0 = -E'0;1, 0,1,0; 

whereas the fourfold degenerate energy level for n = 5 is determined by 

-£■0;!, 1,0,0,0 = -So;i,o,l,0,0 = -E'0;0,1, 0,1,0 = -E'o;0,0,l,l,0- 

However, as can be seen from the figure, for particular values of c it may happen that the 
multiplicity of some eigenvalues is greater than stated in (j3.5p (and hence the number of different 
energy levels is smaller than (I3.4p ) . It is infeasible to obtain analytical expressions for the particular 
values of c for which this will happen. 

Now, we briefly turn our attention to the minimal and maximal energy eigenvalues. For c > 0, 
the values I3k satisfy the following inequalities: 

/3i > /32 > ••• > /3r, fir < Pr+i <■■■< I3n and I3i<l3n, 

and equality between (3r and /3r+i only occurs when n is odd. From this, it immediately follows 
that the maximal eigenvalue of the Hamiltonian in the ladder representation V{p) is given by 
h{l3 + {p — ^)(3n)i and the corresponding eigenvector is wi9, s) with = 1 and Sj = Sj^nijp — !)• It is 
also immediately clear that this eigenvalue is nondegenerate. The minimum eigenvalue is given by 
hpfir- An eigenvector is given by w{9; s) with = and Sj = Sj^rP- When n is even this eigenvalue 
is nondegenerate, but when n is odd it is {p + l)-fold degenerate, since p + 1 is the number of 
compositions of p into two parts. Note that the minimal eigenvalue approaches as c tends to cq. 
This is also illustrated by Figured) 

3.2 Energy eigenvalues in case of fixed wall boundary conditions 

In the case with fixed wall boundary conditions the Hamiltonian ()2.18p becomes, using p. 240 . 

n 

Hpw = eoo + y^/3fcefcfc), 

k=l 

with /3 = ^^1=1 Pk'-, so clearly one has that the basis vectors of the representation are eigenvectors: 

n 

Hfw w{9; s) = h0e + hsk)w{e- s) = hEe,s w{9; s). 

k=l 

The analysis of these eigenvalues turns out to be easier than in the previous case, because the /3j's 
do not satisfy any symmetry relations, but only the inequalities (I2.26p . For c = one recovers the 
same two energy values and degeneracies as before, but in general for c > 0, the energy levels are 
nondegenerate since the f3k do not show any symmetry. Figure [2] illustrates this fact. 

In this case, the maximal eigenvalue of Hpw is fi{(5 + (p — while the minimal eigenvalue is 

given by hpf3n- This follows immediately from the inequalities (12.26p . Of course, these eigenvalues 
are always nondegenerate. 
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(a) n = 4 



(b) n = 5 




Figure 1: (a) The energy levels of the quantum system with periodic boundary conditions for n = 4 
in the representation V{p) with p = 2 and h = uj = I; c ranges from to cq. The vertical axis gives 
the energy values and the numbers next to the levels refer to the multiplicity. When c = there 
are only two energy levels with multiplicities 10 and 4. When < c < cq there are (in general) 
9 energy levels with multiplicity 1, 2 or 3. (b) The same illustration for n = 5. Note that since 
[4/2J = [5/2J the number of different energy levels is (in general) unchanged for < c < cq. 

(a) n = 4 (b) n = 5 




Figure 2: (a) The energy levels of the quantum system with fixed wall boundary conditions for 
n = 4 in the representation V(j)) with p = 2 and h = lo = 1; c ranges from to cq. The vertical 
axis gives the energy values and the numbers next to the levels refer to the multiplicity. When 
c = one recovers the results of the periodic boundary case. When < c < cq however, there are 
(in general) 14 energy levels each with multiplicity 1. (b) The same illustration for n = 5, but now 
there are 20 nondegenerate energy levels. 
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4 On the eigenvalues and eigenvectors of the position operators 



We now turn to the study of the spectrum of the position operators qr-, which are, under the 
solutions (|2.2ip or (j2.24p . quite arbitrary odd elements of the Lie superalgebra gl(l|n). This is 
precisely the topic of [14] where the eigenvalues (and eigenvectors) of the position operators qr 
were determined for arbitrary unitary irreducible representations of the Lie superalgebra 0[(l|n). 
Thus the method developed there will be applied both to the case of periodic and fixed wall 
boundary conditions. 



4.1 Position eigenvalues in the case of periodic boundary conditions 

The position operator g,., given by (j2.7p . can be written as 




h_ 

mn 



n 

(7, e2-^'^/"e,o + 7, e-2-^'^/"eo,) (4.1) 



where we use the notation 

„2 



7j = V^M (i = 1> • • • > ^) and 7 = 7i H h 7„- (4.3) 

The odd g[(l|n) operators 

-.71 n n 

= 4 E 7. e^'^-/"e,o = ^ <e,o, eS = E 4^*^o. (4.4) 
are part of a more general set of operators 



4? = E 4^ and 4) = Y: Uf eoi (1 < J < n), (4.5) 



^^^-Yu^^e^o and <=E^^^'^^* 
/=i 1=1 



where U = {U^^ )i<j<n,i<i<n is a unitary n x n matrix determined by the coefficients in (j4.4p and 



m 



< = ; / , 1 E A...^.^ ^^^- - -^e--0>^)/"e,,,o (4.6) 



(r) (r) 

for J = 1,2, ... ,n — 1. The new operators Ejq and E^j satisfy the same defining relations as 
the elements ejo and eoj; in other words, they generate the Lie superalgebra s[(l|re) [14]. Thanks 
to the introduction of these new operators, one can determine the eigenvalues and eigenvectors of 
the position operators ijr- This follows from the observation that ()4.2p is essentially an element of 
g[(l|l) in the decomposition 1(1 1 n) ^ gl(l|l) ©0l(n — 1). The following result was obtained in [14]: 



Proposition 1 In the representation V{p), all operators qr (r = 1,2, . . . ,n) have the same spec- 
trum. The operator qr has 2p + 1 distinct eigenvalues given by 



The multiplicity of the eigenvalue ±xk ^-^ (" 7^ ) 



±XK = ±\l — {p-K), K = 0,l,...,p. (4.7) 

mn 

n-2+K\ 
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Let us also briefly describe the corresponding orthonormal eigenvectors ipr,±XK,t (with t a mul- 
tiplicity label) corresponding to the eigenvalue ±xk (see [14] for the details). For K ^ p, one 
has: ^ ^ 

1pr,±XK,t = "^^(^'^i' • • • ;*n-l,P - I - K) ± -j=v{0;ti, . . . K), (4.8) 

where ti + ■ ■ ■ + = K. For the eigenvalue {K = p), the eigenvectors read 

1pr,0,t = v{0;ti,. . . ,tn-l,0), h-l ^tn-l=p. (4.9) 

We still need to describe the vectors v{9; t) of V{p) in terms of the basis vectors w{9; s). Essentially, 
the vectors v{9;t) are chosen in such a way that the action of the "new" g[(l|n) elements E^^f} on 
v{9;t) are the same as the action of the "old" gl(l|n) elements ejk on w{9;s). For the new highest 
weight vector, one has: 



„(l;p - 1. 0. . . . , 0) . TJTW^ E(-^)"-"'"'"/C^ 

X 7r'""72 u,p-l-u,0,...,0). (4.10) 
The remaining vectors v{9;t) are given by 

• • • (ii;(J )^-<^-S.-i*^(ii;£))^-*-S}^i*.(i?M)i-^^(l;p - 1,0, . . . ,0); (4.11) 

n~l k 

N = p'-^Uip - -^mt, + l),„^_j;^._^,. (4.12) 

k=l j=l 

with {a)n = a{a + 1) • • • (a + n — 1) the rising factorial symbol. 

Some specific properties of the spectrum of qr will be considered together with those for the 
fixed wall boundary conditions, in the next subsection. The expressions for the eigenvectors will 
be used in the following section, where position probability distributions are studied. 

4.2 Position eigenvalues in the case of fixed wall boundary conditions 

The formal part of the analysis will be rather similar to that of the previous case of periodic 
boundary conditions. But the outcome will show one major difference: the spectrum of qr is now 
dependent (albeit in a simple way) on the position r of the oscillator in the chain. This result is 
completely in accordance with physical intuition as the oscillators in the chain are clearly no longer 
equivalent, since we have two distinguished (and equivalent) oscillators which mark the beginning 
and the end of the chain. 

Under the solution (j2.24p of (j2.19p . the position operator q^., which is given by (j2.14p . becomes 



\ I sin " " ' 



m{n + \)^-^\ n + 



-)7jejo + sin(-^)7jeoj j (4.13) 



y m{n + 1^ 
Here, we used the abbreviations 



{Ki+K:)- (4.14) 



7, = ^andAr,^ = j:sin^(^)7|. (4.15) 
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Comparing with (14. 2p , it is clear that E^^q is obtained from by performing the substitutions 

7j'^7iSin( — — ) (j = l,...,n). (4.16) 

n + 1 

The same substitution can be used for Ejq^ Ej^ and E^^ ^Oj' Then the analysis of the 
eigenvalues and eigenvectors of is determined by the same technique as in the previous subsection. 
In particular: 

Proposition 2 In the representation V{p), the operators Qr (r = 1;2, ...,nj have a spectrum 
depending upon r. The operator Qr has 2p + 1 distinct eigenvalues given by 



±XK = ±J , ^^J p-K), K = 0,l,...,p. (4.17) 
y m[n + 1) 

The multiplicity of the eigenvalue ±xk is ("~^^) • 

The corresponding orthonormal eigenvectors ipr,±XK,t ^'^^ given by the same expressions (j4.8p - 
where one should use the replacement (|4.16p in (|4.1Up (with r = n so that the complex 



exponential reduces to 1) and E^^ £'j^^ in (j4.1ip . 

Although the spectrum of qr is now dependent on the position of the oscillator in the chain 
(i.e. on r) through the constant Ny. in (j4.17p . one still has the expected symmetry that the spectra 
of qr and ^n-r+i coincide, since = A^^_j._|_i. 

The dependence on r of the spectrum of qr is completely determined by Nr, which on its turn 
depends on the coupling constant c. In Figure O we will plot some of these eigenvalues as a function 
of c. For these plots, we use m = h = uj = 1, and choose K = p — 1 (then the eigenvalue expression 
is independent of p). We plot the cases of periodic boundary conditions and of fixed wall boundary 

conditions in two different figures. Thus for periodic boundary conditions, we plot the value 



as a function of c; for fixed wall boundary conditions, we plot the values -i/ — rr-^r as a function of 



C. 

This figure, together with other numerical experiments, suggests the following about the range 
of the spectrum, independent of any measurement probabilities: 

• When the coupling constant c is fixed one has that A^i < N2 < • • • < -^[2.] . This means 
that the spectrum of oscillators close to the wall is more centered around their equilibrium 
positions than for oscillators in the middle of the chain. This is an intuitively clear result as 
oscillators in the middle of the chain do not "feel" the walls as much as oscillators close to 
the walls do. 

• Viewed as a function of c, A'^i is decreasing. When the coupling constant c increases, it is 
clear that the movements of the first oscillator will become more restricted, resulting in a 
spectrum closer to its equilibrium position. 

• The spectra of position operators associated with oscillators in the middle of the chain seem 
to become very similar to one another (as n increases). This is in line with intuition as these 
oscillators are more or less equivalent with respect to their distance to the wall. Related 
to this fact one sees that the spectrum of an arbitrary position operator associated with an 
oscillator from a chain with periodic boundary conditions is also very similar to the spectrum 
of an oscillator in the middle of the chain with fixed wall boundary conditions. 
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Figure 3: Indication of the spectrum of various position operators for the two chains, as a function of 
the couphng constant c and with m = H = u = 1. The horizontal axis in these figures represents c. 
In the left column n = 4, while in the right column n = 6. In the top row ^y^jjn, i.e. the eigenvalue 
for an arbitray position operator (r = 1, . . . ,n) in a chain with periodic boundary conditions with 

K 



p — 1, IS snown. i ne ngures m tne secona row snow y ^^A^^j with r indicated next to the 
graph. In other words, they give the eigenvalue of the position operators in a chain with fixed wall 
boundary conditions (again with K = p — 1). 
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• For oscillators other than the first and the last, at first the spectrum moves closer to the 
equilibrium positions, but then seems to widen again. This effect could be related to an 
increasing effect of "collective" motions (lattice vibrations) as the coupling constant becomes 
larger. 

In the above figures, we have concentrated on a particular K-value, i.e. on a particular eigen- 
value. Of course, the complete spectrum of each operator is easy to describe: it is simply a 
constant times ±\/p — K, K = 0, 1, . . . This is a simple distribution that will be depicted in the 
figures of the following section. 



5 Position probability distributions in the ladder representations 

In this section, we shall establish some facts about the position probability distributions of the 
two systems. It is a well known fact (postulate) of quantum mechanics that when measuring an 
observable, the measurement always yields an eigenvalue of the (self-adjoint) operator associated 
with that observable. The probability of measuring a certain eigenvalue when the system is in a 
certain state is determined by the expansion of that state in terms of (orthonormal) eigenvectors 
of the operator at hand. 

We now assume that our system is in a stationary state w{9] s) and we wish to determine the 
probabilities of obtaining the different eigenvalues xk of an operator g.^ . In the ladder representation 
V{p), the operator has 2p + 1 distinct eigenvalues ^xk = \{p — K), where Q < K < p 
and with a constant that in the case of periodic boundary conditions is independent of r, but 
does depend on r in the case of fixed wall boundary conditions. The multiplicity of the eigenvalue 
zbxi^ is ("'~^^)) and the orthonormal eigenvectors ipr,±XK,t (with t a multiplicity label) have been 
given in the previous section. Each eigenvector 'tpr,±XK,t can be expanded in terms of the stationary 
states w{6; s): 

A,±XK,t = Yl ^rl±XK,t^^^'^ 

e,s 

Using (for the case of periodic boundary conditions) equations (|4.1U|) and (j4.11|) the coefficients 

s 

C^'_l_j.^ f in this expansion can be computed explicitly (the same holds of course for the case of 
fixed wall boundary conditions). Using orthonormality of the eigenvectors and the basis vectors, 
we immediately have that 

When the quantum system is in the fixed stationary state w{9^ s), the probability of measuring for 
qr the eigenvalue ±xk is given by: 

P{e,s,r,±XK)= Y |<±wf- (5-1) 

ti + --+t„-i=K 

From ()4.8p it is immediately clear that 

P{e,s,r,XK) = P{e,s,r,-XK). (5.2) 

When p = 1, one can determine explicit expressions for the vectors v{(j);t). In this case when 
^ = and there is exactly one tj = 1, we introduce the shorthand notation v{0; V) for this vector. 
In the case of periodic boundary conditions, the position probabilities do not depend on the position 
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r of the oscillator in the chain, so we work with r = n so that all complex exponentials occurring 
in the various expressions reduce to 1. In this case, the highest weight vector is given by 

z;(l;0,...,0) =u;(l;0,...,0), 
and using induction on k, (since v{0; 1^^^) = E^-^ k'^i^'i l'^)) '^^^ prove that 

k 



v{0; 1'=) 
^(0; 1") 



1 



7i7fc+i 



7? + • • • + ll+i 'j=i + • • • + 7i 



(0; P) -Jjf + ...+ 71 y;(0; i^'+i)), when k + 



n 



^7,u;(0;l^). 



We can now immediately give the expansion of the eigenvectors of g„ in terms of the basis vectors 
w{d;s): 



n-l 



'4'n,±xo,t 



-^wil; 0, . . . , 0) ± -1= ^7^.u;(0; 1^), t = (0^^). 



7? + ••• + 7^+1 'j=i + + 



^0;P)- J7? + ---+7i«^(0;l^+i)), t = 



with \ < k < n — 1. From this and (jS.ip it immediately follows that 



P{9,s,r, xo) 



and hence 



1^ when0 = O,Sfc = l 
^ when ^ = 1, 



P{e,s,r,x,) = { 7 when^ = 0,.. = l 
when = 1. 



Also for p = 2 we have been able to determine the position probabilities. They are as follows: 



2^ 



P{e,s,r,xo) = < 



when 9 = 0, Sk = 2 



P{e,s,r,xi) = < 



when 9 = 0, Sk = si = 1 

1 V 

2^ when 6* = = 1, 



Till^ when e = 0, Sfc = 2 

1 (7-7|-7?)(7g+7?)+(7F7?)^ 



7 



1 7-7fc 

2 7 



when 6* = 0, Sfc = = 1 
when 6* = 1, Sfc = 1, 



and 



P{9,s,r,X2) = < 



when 9 = 0, Sk = 2 



(7-7F7?)7+27i7F ^^^^ ^ ^ ^ ^ ^ 





when = 1, Sfc = 1. 



In the case of fixed wall boundary conditions, these position probabilities are given by the same 
expressions, subject to the substitutions (|4.16p . So, as was to be expected, the position probabilities 
do depend on the position of the oscillator in the chain. 
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For p = 1 and p = 2 it was possible to compute these position probabilities analytically, but for 

S 

p > 2 this becomes infeasible. On the other hand, since all coefficients C^'-j-^, j are known explicitly 
(in all cases and for any p), we can numerically compute all position probabilities. We will now 
examine the plots of some of these position probability distributions. 

Let us first consider the case of periodic boundary conditions. The position probabilities are 
independent of r, as the system is completely symmetric, so let us take r = 1. We will plot the 
values 

P{d,s,l,±XK), K = 0,l,---,P 

for certain values of 9 and s. In other words we plot the position probability distribution function 
when the system is in a fixed stationary state ^^{O; s). Let us consider an explicit example, say n = 4 
(four coupled oscillators) and p = W (so each position operator has 21 distinct eigenvalues). We 
will plot the position probability distributions for the ground state (this is the state w{0; 0,p, 0, 0) = 
w{0; 0, 10, 0, 0)) and for the most excited state (this is the state 'w{l; 0, 0, 0,p — 1) = w{l; 0, 0, 0, 9)). 
These distributions are given in Figure [H for some c- values. 

Let us make a number of observations on these distributions. When the system is in the ground 
state, the probability distribution function of each position operator is symmetric around its equi- 
librium position and unimodal. Of course it is also discrete (as we are working in finite-dimensional 
representations). As the coupling constant c increases, the "peak" around the equilibrium position 
is sharper. In other words, a)s the coupling constant becomes larger, the oscillators are more likely 
to be close to their equilibrium position when the system is in its ground state. 

When the system is in its most excited state, the position probabilities are quite different. The 
probability of finding the oscillator in its equilibrium position is zero. On the other hand, there are 
certain peaks away from the equilibrium position. As c increases, these peaks are further away from 
the equilibrium position. In other words, as the coupling constant becomes larger, the oscillators 
are more likely to be further away from their equilibrium position when the system is in its most 
excited state. 

Note that in this figure one also observes the fact that the range of the spectrum of the position 
operators is dependent on the coupling constant c. When talking about probabilities of being 
further or closer to the equilibrium position we regard this relative to the discrete spectrum of 
2p + 1 values (the middle one being the equilibrium position). 

Let us now consider the case of fixed wall boundary conditions. The situation is rather different, 
as the position probabilities are depending on r. To see the r-dependence, we will plot position 
probability distribution functions for r = 1 (the oscillator just next to the fixed wall) and for r = 3 
(an oscillator away from the wall). Again, we will plot the values 

P{9,s,r,±XK), K = 0,l,...,p 

for certain values of 6 and s, i.e. when the system is in a fixed stationary state w{9;s). As an 
explicit example, take n = 6 (six coupled oscillators) and p = 7 (so each position operator has 15 
distinct eigenvalues). We will again plot the position probability distributions for the ground state 
(this is the state t(;(0; 0, 0, 0, 0, 0,p) = t(;(0; 0, 0, 0, 0, 0, 7)) and for the most excited state (this is 
the state w{l;p — 1, 0, 0, 0, 0, 0) = wi^l; 6, 0, 0, 0, 0, 0)). These distributions are given in Figure [5] for 
r = 1 and in Figure [6] for r = 3, for a number of c- values. 

Let us again make a number of observations on these distributions. When the system is in the 
ground state, the probability distribution function of each position operator is symmetric around its 
equilibrium position and unimodal. For fixed c > 0, if the oscillator is closer to the wall (r = 1) the 
peak of the distribution function around the equilibrium position is sharper than for an oscillator 
further away from the wall (r = 3). In other words, the oscillators close to the wall are closer 
to their equilibrium position than those further away from the wall. As the coupling constant c 
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Figure 4: Position probability distribution function for the position operator qi, in the periodic 
boundary case, when n = 4 and p = 10. In the three rows, c = 0.1, c = 0.4, c = 0.8. The plotted 
value is P{9, s, 1, ±xk) for each of the 21 eigenvalues ±xk {K = 0,1, . . . , 10) of qi. This is given 
for the case when the system is in the stationary state w{9;s) corresponding to the ground state 
(minimum energy) in the left column and in the right column when it is in the stationary state 
w{9;s) corresponding to the most excited state (maximum energy). 
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Figure 5: Position probability distribution function for the position operator qi, in the fixed wall 
boundary case, when n = 6 and p = 7. In the three rows, c = 0.05, c = 0.2, and c = 0.4. The 
plotted value is P{9, s, 1, ±xk) for each of the 15 eigenvalues ±xk [K = 0, 1, . . . , 7) of gi (next 
to the wall). In the left column, this is given for the case when the system is in the stationary 
state w{9; s) corresponding to the ground state (minimum energy) and in the right column for the 
case when the system is in the stationary state w{9] s) corresponding to the most excited state 
(maximum energy). 
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Figure 6: Position probability distribution function for the position operator q^, in the fixed wall 
boundary case, when n = 6 and p = 7. In the three rows, c = 0.05, c = 0.2, and c = 0.4. The 
plotted value is P{9, s, 1, ±xk) for each of the 15 eigenvalues ±xk {K = 0, 1, . . . , 7) of gi (next 
to the wall). In the left column, this is given for the case when the system is in the stationary 
state w{9; s) corresponding to the ground state (minimum energy) and in the right column for the 
case when the system is in the stationary state wiO^s) corresponding to the most excited state 
(maximum energy). 
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increases, the "peak" around the equihbrium position also becomes sharper, both for r = 1 and 
r = 3. In other words, as the couphng constant becomes larger, the oscillators are more likely to 
be close to their equilibrium position. 

When the system is in its most excited state, the position probabilities are rather different. 
The probability of finding the oscillator in its equilibrium position is zero. For c > 0, there are 
certain peaks away from the equilibrium position, both for r = 1 and for r = 3. Close to the wall 
(r = 1), these peaks are closer to zero than away from the wall (r = 3). So also in this most excited 
state, the oscillators close to the wall "oscillate less heavily" than those away from the wall. As 
c increases, these peaks are further away from the equilibrium position, both for r = 1 and for 
r = 3. So the oscillators are more likely to be further away from their equilibrium position when 
the coupling constant increases. 

6 Conclusions 

We have examined properties of noncanonical solutions of two quantum systems: chains of coupled 
harmonic oscillators with periodic boundary conditions or with fixed wall boundary conditions. 
These new solutions arise from an approach as a WQS, allowing more classes of solutions than just 
the canonical one. 

For the solutions examined here, the position and momentum operators are (odd) elements of the 
Lie superalgebra 0l(l|n). The physical properties of the system then follow from the representations 
of gl(l|n) considered. Here, we have introduced the rather simple ladder representations V(j>), a 
class of unitary irreducible representations of fl[(l|n) (or rather, of its compact form u(l|n)). For 
these representations, we have determined the energy spectrum, which was rather easy due to the 
simple action of the Hamiltonian operator in the standard basis of V{p). We have also determined 
the spectrum of the position operators; this task was more difficult because of the more complicated 
action of these operators in the basis of V{p). The techniques developed in [14] allow to construct 
explicitly the eigenvectors of the position operators, for both systems under consideration. 

The analysis of the spectrum of the position operators, and their probabilities when the system 
is in a certain stationary state, lead to interesting properties. The spectrum is discrete, centered 
around the equilibrium position; the number of possible position values depends on p (it is 2p+ 1). 
The width of this discrete support depends on the coupling constant c, and - only in the case 
of fixed wall boundary conditions - on the order r of the oscillator in the chain of n oscillators. 
The position probability distributions, discussed in detail in the previous section, have properties 
similar to those of a classical system of coupled oscillators. 

The current paper was still dealing with solutions arising from the Lie superalgebra £|[(l|n), 
although a different class of representations was considered than in [7] . As indicated in a previous 
paper [7], also other types of solutions for the compatibility conditions exist, for example in terms 
of the Lie superalgebra osp(l|2n). It would be interesting to investigate other solutions in terms of 
this orthosymplectic Lie superalgebra, even though the analysis is expected to be rather difficult. 
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